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We evaluate the prospects for finding evidence of dark matter production at the Large Hadron 
Collider. We consider WIMPs and superWIMPs, weakly- and superweakly-interacting massive parti¬ 
cles, and characterize their properties through model-independent parameterizations. The observed 
relic density then implies lower bounds on dark matter production rates as functions of a few pa¬ 
rameters. For WIMPs, the resulting signal is indistinguishable from background. For superWIMPs, 
however, this analysis implies significant production of metastable charged particles. For natural 
parameters, these rates may far exceed Drell-Yan cross sections and yield spectacular signals. 
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The energy density of non-baryonic dark matter in the 
Universe is now known to be 0 

= 0.112 ±0.009 , (1) 

where Udm is this energy density in units of the critical 
density, and h ~ 0.71 is the normalized Hubble param¬ 
eter. With accompanying constraints, this implies that 
non-baryonic dark matter makes up about a quarter of 
the energy density of the Universe. Its microscopic iden¬ 
tity is at present unknown, however, and is one of the 
outstanding questions in basic science today. 

An exciting possibility is that particles that make up 
some or all of dark matter may be produced at high- 
energy colliders, such as the Tevatron at Fermilab or, 
beginning in 2007, the Large Hadron Collider (LHC) at 
CERN. Such prospects are particularly promising if dark 
matter is composed of WIMPs or superWIMPs, weakly- 
or superweakly-interacting massive particles, since these 
scenarios require new particles with masses near M^eak ~ 
100 GeV, the scale to be probed in detail at the LHC. 
Of course, the discovery of dark matter signals requires 
not only that dark matter be produced at colliders, but 
also that it be produced with rates and signatures that 
allow it to be distinguished from background. Remark¬ 
ably, however, WIMP and superWIMP production rates 
are also constrained by cosmology, because dark mat¬ 
ter densities are determined (in part) by thermal freeze 
out in these scenarios. Since WIMPs and superWIMPs 
cannot have relic densities in excess of they must 

have annihilated efficiently in the early Universe and con¬ 
sequently must be produced efficiently at colliders. 

In this paper, we analyze this argument quantitatively. 
We begin by characterizing the properties of WIMPs and 
superWIMPs through model-independent parameteriza¬ 
tions. Our analysis is sufficiently general to accommo¬ 
date concrete realizations of WIMP or superWIMP dark 
matter in supersymmetric, extra-dimensional, and many 
other frameworks. Using the observed relic density, we 
then determine, as functions of the few parameters en¬ 


tering the analysis, the minimal dark matter production 
cross section and its observable consequences. This ap¬ 
proach was applied in Ref. 0 to obtain conservative esti¬ 
mates for WIMP cross sections at the International Lin¬ 
ear Collider (ILC). The results implied challenging rates 
for discovery. In this study, we focus on the more pressing 
case of the LHC and also analyze superWIMP scenarios, 
for which our conclusions will be far more promising. 

In both the WIMP and superWIMP dark matter sce¬ 
narios, the process of freeze out plays a large role in de¬ 
termining the dark matter relic density. In WIMP sce¬ 
narios, the dark matter is a stable, neutral WIMP with 
mass mwiMP A^veak and pair annihilation cross section 
a A '^weak-^weak- Beginning in thermal equilibrium in 
the early Universe, WIMPs annihilate until they become 
too dilute to find each other and “freeze out” at tem¬ 
perature Tp. Given the expansion rate determined by 
Mpi ~ 1.2 X 10^® GeV, along with mwiMP and cta given 
above, the thermal relic density Uwimp^^ is automati¬ 
cally near the observed of Eq. (P). In contrast to 

other dark matter scenarios, there is no need to introduce 
new energy scales to obtain the desired relic density. 

SuperWIMP scenarios also begin with a weakly- 
interacting particle, which we denote L, with mass rrip 
Mweak and cta ^ “Lak-^welk freezes out with den¬ 
sity Ylph'^ ^ In contrast to WIMP scenarios, 

however, L particles are not stable, but metastable, and 
they ultimately decay to superWIMPs, which form the 
dark matter we observe today 0 ■ SuperWIMPs are neu¬ 
tral and stable, but their interactions are much weaker 
than weak, typically gravitational. The resulting dark 
matter density is, then, 

SiswiMptf = , (2) 

rriL 

If, as is often natural, mswiMP mp, superWIMPs 
are also produced with relic densities of the right or¬ 
der of magnitude. A priori, L particles may be ei¬ 
ther electrically charged or neutral. However, in the 
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well-motivated cases of gravitino and Kaluza-Klein gravi¬ 
ton superWIMPs i II i, the neutral case is typi¬ 
cally excluded by contraints from Big Bang nucleosyn¬ 
thesis ^ 1^ , and the most motivated L particles are 
charged sleptons and KK leptons. In this study, we there¬ 
fore consider only L particles with charge ±1. 

To determine the production rates at colliders, we must 
first find what annihilation cross sections are implied by 
the observed relic density Let X denote a generic 

particle that freezes out, either the WIMP in WIMP sce¬ 
narios, or L in superWIMP scenarios. The total XX 
annihilation cross section is 


o-tot = s) . (3) 


where •\/i is the center-of-mass energy, and i,j are par- 
tons. We parameterize the thermal average of Utot as 


(atotvx) = aant'i" + 0(4”^') = ^ox-^ + 

(4) 

where vx is the relative velocity of the initial X and X 
particles in their center-of-mass frame, and x = mx/T = 
Gjvx- This expansion is valid in the common case where 
annihilation is dominated by a single angular momentum 
component (S'-wave for n = 0, P-wave for n = 1, etc.). 
It is necessarily valid at freeze out, where a; ^ 25 and 
Vx ^ but, of course, breaks down if vx is near 2. 

With this parameterization, the annihilation cross sec¬ 
tion is related to the relic density through 0 


^xh^ ^ 1.07 X 10^ GeV~^ ^ 

v^Mpi ao 

1 

° — 1V 45 g mxMpi 


(5) 

( 6 ) 


where xp = mx/Tp^ is the effective number of 
relativistic degrees of freedom at freeze out, g is the 
number of X degrees of freedom, and c is defined by 
Y(xp) = cYpq(xp)^ with Y(x) the X number density 
per entropy density and lEQ(a;) its value if X had re¬ 
mained in thermal equilibrium. We set — 1 = n -|- 1, 
which reproduces numerical results to within 5% 0 . 

Given the annihilation cross section for X particles, 
their production rate at colliders is fixed by the principle 
of detailed balance, assuming time reversal symmetry 0. 
Neglecting the i,j parton masses, we find 


a{ij XX; s) 


4(25, + 1)(25, + 1) ^ 


*j; s) 


giji‘2Sx + 1 )^ 
4(25, -h l)(25j -h 1) 




2n+l 

X 


,(7) 


where Kij = a{XX ij; s)/crtot, Vij is ^ if f and j are 
identical and 1 otherwise, 5 denotes spin, and all cross 
sections include averaging and summing over initial and 
final state spins, respectively, but do not include color 


factors. Note that Eq. 0 has been used, and so the fi¬ 
nal expression of Eq. 0 is not trustworthy for vx near 
2. Equations 0-0 determine the minimum production 
cross section, given the observed relic density, as a func¬ 
tion of a few parameters that characterize the properties 
of X: its mass mx and spin Sx, its dominant annihila¬ 
tion channel n, and the dynamical parameters k^-. 

We now turn to the case of superWIMP dark mat¬ 
ter produced in L decays. The L lifetime is extremely 
long (for gravitational decays, it is of the order of hours 
to months), and so L particles appear as stable charged 
particles in colliders. Eor the LHC, with i/i = 14 TeV, 
the minimum cross section for L pair production is 

1 

a{pp^ L+L-;s)= du ^ 

X Ey [fl ix)fl (u/x) + fl {x)fl (u/x)] 
xaigidj L^L~;us) , (8) 

where ff are proton parton distribution functions, and 
d{q^qj L+L-;us) = ^ cr{q^qj L+L~;us) (9) 

^ color 

is the color-averaged parton-level cross section, where 
Nc = 3, and the right-hand side is determined by Eq. 0 
with X replaced by L. The upper limit of integration for 
u forces vp < Umax- We choose = 2, so that the 
parametrization of Eq. 0 may reasonably be expected 
to be valid in the region of integration, and we conserva¬ 
tively neglect all contributions from vp > Umax- We also 
neglect subleading contributions from (loop-induced) gg 
fusion and three-body processes qg —>■ L~q. Tevatron 
cross sections are determined by replacing one proton 
with an anti-proton and setting = 2 TeV. 

To distinguish the metastable L signal from back¬ 
ground, we require that both L+ and L~ have pseudo¬ 
rapidity \r]\ < 2.5 and velocity f3 < 0.7, so that both 
tracks will be detected with ionization —dE/dx more 
than double minimum-ionizing. The j3 and vx require¬ 
ments are correlated but independent, since (3 is in the 
lab frame, and vx is in the parton center-of-mass frame. 
These cuts, together with the requirement of isolated 
tracks in events free of hadronic activity, should leave 
the signal essentially background-free. The event rates 
depend weakly on the i] cut, dropping by about 20 % 
when requiring I 77 I < 0.5. For /3 < 0.6, the event rate 
drops by a factor of 2 to 5, depending on mp. 

In Fig. ^ we show the minimum cross sections for L 
pair production at the Tevatron and the LHC as func¬ 
tions of mp for both 5- and P-wave annihilation, assum¬ 
ing scalar L particles, velocity-independent Kqq =0.2 for 
q = d,u,s,c,b, and mswiMP= 0.6. As can be seen, 
P-wave annihilation in the early Universe is suppressed 
by u| relative to 5-wave, and so must be compensated 
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FIG. 1: Minimum Tevatron and LHC cross sections derived 
from cosmology, along with the Drell-Yan cross section, for 
L^L~ production in the superWIMP scenario as functions of 
niL- Each event has two highly-ionizing tracks with velocity 
(3 < 0.7 and psuedorapidity \ri\ < 2.5. We assume Sl ~ 0, 
mswiMp/m-i = 0.6, and Kqq = 0.2 for q = d,u, s, c, b. 

by larger (7an, leading to larger minimum collider rates. 
D- and higher wave annihilation will imply even larger 
minimum rates. The dependence on the other parame¬ 
ter assumptions is that the cross sections scale linearly 
with {2Sl + and, to an excellent approximation, 

toswimp/wl. The L particles are assumed to be pro¬ 
duced isotropically in the parton center-of-mass frame. 
We have checked that the results are insensitive to this 
assumption, varying by less than 10% for alternative dis¬ 
tributions, such as sin^ 9 and (1 ± cos0)^. 

The cosmologically constrained L production cross sec¬ 
tions shown in Fig. ^ are significant. For the LHC, 
even for tol ~ 1 TeV, cross sections as large as 0.1 fb 
are predicted for P-wave annihilation. For comparison, 
the Drell-Yan production cross section a{pp —>■ "f, Z 
L+L“) is also shown in Fig. Q] At the LHC, for the pa¬ 
rameters chosen, the minimum cross sections we have 
derived typically exceed the Drell-Yan cross sections; 
for P-wave annihilation, they are bigger by factors of 
2 to 50, depending on mz,. Note that these minimum 
cross sections are parameter-dependent; no absolute min¬ 
imum can be derived, as, for example, these cross sec¬ 
tions may be made arbitrarily small by taking mswiMP 
to zero. However, if the fact that HswiMp/i^ is around 
is not merely a coincidence, it is natural to expect 
wswiMP ^ "rriL, and it is significant that the predicted 
cross sections may be large for such natural choices. 

The analysis above may be adapted easily to the case of 
the ILC, where the role of Kqq is played by Kg+e- ■ For the 
ILC, /3 and nmax/2 are nearly identical; the derivation of 
collider cross sections from the relic density is therefore 
reliable only when L particles are produced with (3 < 
0.7. As examples, consider Sl = 0, mswiMp/»7iL = 0.6, 
Kg+g- = 0.1, and L particles produced with (3 = 0.7. 
For both ^/s = 500 GeV and 1 TeV, we find minimum 



FIG. 2: Discovery limits for dark matter events in the super¬ 
WIMP scenario at the Tevatron with L = 30 fb“^, LHC with 
L = 1 ab“^, and ILC with L = 1 ab~^. We require 10 events 
with two doubly minimum-ionizing tracks for discovery. For 
the ILC, the beam energy is assumed to be ^/s ~ 2.8 tul 
so that L particles are produced with /3 = 0.7. We assume 
Sl = 0, P-wave annihilation, and mswiMp/mn = 0.6. The 
reach in k scales linearly with (25'^-I-1)“^ and, to an excellent 
approximation, (mswiMp/m-i,)”^. 

cross sections of roughly 4 fb and 60 fb for S'-wave and 
P-wave annihilation, respectively. These imply hundreds 
to thousands of background-free events per year. 

In Fig. [3 we give discovery limits in the (mz,, k) plane, 
where k = Kqq for the Tevatron and LHC, and k = Kg-t-g- 
for the ILC. Given the cuts described above, we ex¬ 
pect the signal to be background-free, and so require 
10 signal events for discovery. Even assuming a highly 
optimistic luminosity, the Tevatron can see the mini¬ 
mum dark matter signal only for itil ^ 100 GeV and 
near maximal Kqq. The LHC does much better, probing 
Kqq > 3 X 10“^ for rriL ^ 100 GeV, and mz, < 1.2 TeV 
for Kqq ~ 0.2. Finally, the ILC will provide phenomenal 
coverage down to Kg+g- ~ 10“'* for all m^. The process 
L+P“ ^ e+e" may be absent at tree-level; in fact, all 
processes L^L~ —> // may be suppressed if annihila¬ 
tion to Higgs bosons dominates. However, even in these 
cases, since L must be weakly coupled to some standard 
model particles to explain flLh'^ ~ we expect 

Kg-i-g- ~ 10“^ to be generated at loop-level. Of course, 
this requires kinematically accessible L pairs. This is rea¬ 
sonable for the lower range of itil plotted, but requires 
later stages of the ILC program for the upper range. 

We now consider WIMP dark matter scenarios. WIMP 
pair production is invisible at colliders. Here we con¬ 
sider the monojet signal of WIMP pairs produced in as¬ 
sociation with a gluon or quark. The cross section for 
pp XXj is not simply related to the cross section for 
pp XX. However for collinear or soft jets, these rates 
are related by splitting functions . The color-averaged 
parton-level differential cross sections are, then, 
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TABLE I: The minimum monojet signal S from pp —> XXj 
at the LHC in the WIMP dark matter scenario. We assume 
scalar WIMPs with mass 100 GeV and P-wave annihilation 
and require jets to have sin Slab >0.1 and as indicated. 
Also given are the standard model background B and the 
significance S/y/B, assuming integrated luminosity 1 ab“^. 

5 _ B _ s/Vb 

30 GeV 18.6 fb 1300 pb 0.51 

100 GeV 4.1 fb 130 pb 0.36 


Fq^g{z, 0) a{qq XX- (1 - z)s) 


( 10 ) 


dz dcos0 


m) XXq-, a) 

d{qq^ XX-, (1 - z)s) , (11) 


25 , + 1 . 


25„ + l 


where the splitting functions are [9j 

4 Os 1 + (1 - z)^ 1 




N, 


sin 0 


( 12 ) 


P (~ n - ^ asz^ + {l-zf 

— 1 TT sin 0 


with identical expressions for q q. Here i ^ jk means 
that initial state parton i radiates parton fc, which be¬ 
comes the final state jet. The parameter z = Ek/Ei 
varies from 0 to 1 — 4m^/[s(l — u^ax/^)] and 0 is the 
angle between i and fc; both are defined in the parton 
center-of-mass frame. Given these parton-level results, 
the LHC color-averaged differential cross section is 


f 

TT^ 


d{pp XXg, XXq, XXq-, a) 
du 


dz dcos0 

4m^ 


^ dr ^ 




dzd cos 0 


a(i(j ik) XXk]ua) . 


(14) 


To ensure that the monojet events are detectable, we 
require the jets to have sin^iab > 0.1 and pr > p™™. 
The factorization of Eqs. m and holds formally 
only in the limit of soft or collinear jets, but it has been 
shown to be reasonably accurate even away from these 
limits in the region we have included |^. In contrast to 
the super WIMP case, where the background is negligible, 
monojet events in the WIMP scenario suffer from a huge 
irreducible background from pp —>■ vkj. At the ILC, the 
analogous single photon signal may be improved by an 
additional cut on the photon energy |^. Unfortunately, 
this approach is not effective at the LHC because the par- 
ton center-of-mass energy is not fixed. Table Ogives cross 
sections for the monojet signal for mwiMP = 100 GeV 
and two values of p™'". The background, with the identi¬ 
cal cuts implemented using the simulation package COM- 
PHEP [ 13 , is also given. Although thousands of WIMP 


monojet events are expected given integrated luminosity 
1 ab“^, the overwhelming background leads to very small 
S/'/B, making discovery extremely difficult. 


In summary, if stable WIMPs or super WIMPs exist, 
requiring that they not overdose the Universe implies 
efficient dark matter production rates at colliders. Us¬ 
ing model-independent parametrizations, we have deter¬ 
mined lower bounds on dark matter production rates as 
functions of a few parameters characterizing WIMP and 
superWIMP properties. For WIMP dark matter, the 
XXj signal is swamped by background. On the other 
hand, for natural parameters in the superWIMP scenario, 
the derived rate for the production of two metastable 
charged particles may be much larger than the Drell-Yan 
cross section and implies spectacular signals at the LHC 
and, if kinematically accessible, the ILC. These super¬ 
WIMP results imply promising prospects not only for 
detection of dark matter signals, but also for detailed 
studies 0 of dark matter at future colliders. 
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